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ABSTRACT 

The Swift Burst Alert Telescope (BAT) is discovering interesting new objects 
while monitoring the sky in the 14-195 keV band. Here we present the X-ray prop- 
erties and spectral energy distributions for two unusual AGN sources. Both NVSS 
193013-F 341047 and IRAS 05218-1212 are absorbed, Compton-thin, but heavily ob- 
scured (Nh ~ lO^^cm"^), X-ray sources at redshifts < 0.1. The spectral energy dis- 
tributions reveal these galaxies to be very red, with high extinction in the optical and 
UV. A similar SED is seen for the extremely red objects (EROs) detected in the higher 
redshift universe. This suggests that these unusual BAT-detected sources are a low- 
redshift {z « 1) analog to EROs, which recent evidence suggests are a class of the 
elusive type II quasars. Studying the multi-wavelength properties of these sources may 
reveal the properties of their high redshift counterparts. 

Subject headings: surveys, X-rays: galaxies, galaxies: active 



1. Introduction 

Super massive black holes are believed to exist at the centers of most massive galaxies. While 
a connection between the central black hole and larger host galaxy has been observationally proven 
(the M - a relation), there are still many unanswered questions as to how the central black hole 
affects the host galaxy's formation and evolution. Since previous studies of active galactic nuclei 
(AGN) have relied heavily on optically or X-ray sele cted samples, they missed a significant fractio n 



(ss 20%) of AGNs - the heavily obscured sources (|Ueda et al.l l2007l : IWinter et al.l l2008l . I2009al 1. 
Thus, our understanding of the AGNs and their connection to the host is missing a potentially key 
parameter. 



^Department of Astrophysical & Planetary Sciences, University of Colorado, UCB 391, Boulder, CO 80309, USA 
^Department of Astronomy, University of Maryland, College Park, MD 20742, USA 
^Hubble Fellow 



- 2 - 



Surveys of active galactic nuclei are typically dominated by two selection effects: dilution by 
starlight from the host galaxy and obscuration by du st and gas in the host galaxy and the AGN 
itself (see iHewett fc: Folta Il994l and iMushotzkvl |2004| ) . These factors previously kept an unbiased 
AGN sample from reach. However, with the capabilities of Swift's Burst Alert Telescope (BAT), 
this has changed. The BAT surveys the sky in the hard X-ray range of 14 - 195 ke V. In the first 
22 m onths of the survey, BAT detected ~250 AGN with a median redshift of 0.03 (iTueller et al 



2010l ). These AGN were selected purely by their hard X-ray flux, and thus, are not affected by 
obscuration by gas and dust below ~ lO^^cm"^, which prevents them from being easily detected 
in optical or soft band X-ray surveys. 

In the BAT sample there is no bias against detecting heavily obscure d but still Comp ton thin 
AGNs, thus the sample offers a unique test of the unified model of AGNs (j Antonuccil 1 19931 ) . In the 
unified model, differences between broad (type 1) and narrow- line (type 2) AGNs a re attributed 
to differences in the viewing angle to the central, obscured region (|Osterbrocklll978l ). Therefore, 
there should be no differences between the X-ray and optical AGN classifications. 

However, the X-ray/optical properties of two sources out of the ~ 250 detected do not fit the 
most basic prediction of the unified model. Both NVSS 193013-^341047 and IRAS 05218-1212 show 
broad optical lines indicative of a type 1 - 1.5 optical source, but the low resolution Swift XRT 
spectra indicated that the sources were likely highly absorbed. Further, the combined BAT and 
XRT X-ray spectra indicated that the sources were mildly "Compton-thick" , with column densities 
of Nh > 1.4 X lO^'^cm-^. The spectral energy distributions (SEDs) of these sources are extremely 
red, showing very heavy extinction levels in the UV and soft X-rays (< 2keV). Such SEDs have 
only previously been seen for the high redshift {z > 1) extremely red objects (EROs). In order to 
determine the X-ray properties of these sources, we obtained X-ray follow-ups with XMM-Newton. 
In this paper, we present both the X-ray and broad-band properties of these two unusual sources 
detected in the Swift BAT survey. 



2. X-ray Spectra 

2.1. Observation Details 

XMM-Newton follow-ups were acquired in the 0.1-10 keV band for both NVSS 193013-^341047 
and IRAS 05218-1212. Basic details of the sources and the XMM-Newton observations are included 
in Table [H The spectra were processed and extracted using the Science Analysis System (SAS) 
version 9.0. We followed the process for reduction outlined in the XMM-Newton ABC Guid^. First, 
we created calibrated photon event files for the EPIC MOS and pn cameras using the observation 
data files (ODFs) with the commands emchain and epchain, respectively. 



^The ABC Guide is available online at |http : / /heasarc . nasa . gov/ docs/xmm/ abc/ 1 
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The events tables were filtered using standard filters as outlined in the ABC Guide. For 
MOS data, good events constitute those with a pulse height in the range of 0.2-12 keV and event 
patterns that are characterized as 0-12 (single, double, triple, and quadruple pixel events). For the 
pn camera, only patterns of 0-4 (single and double pixel events) are kept, with the energy range 
for the pulse height set between 0.2 and 15 keV. Bad pixels and events too close to the edges of 
the CCD chips were rejected using the most stringent selection expression, FLAG == 0, for both 
data sets. Based on an inspection of the light curve, we found no evidence of flaring in either 
observation. Therefore, we did not apply any temporal filtering. 

The spectra of the sources were extracted using the SAS task evselect. The spectra were 
extracted in a circular aperture centered on the source. Additionally, background spectra were 
extracted with a circular aperture of the same size, in a region nearby the source on the same 
CCD chip, but free of additional X-ray sources. Response files were generated using the SAS tasks 
rmfgen and arf gen. The spectra were then binned by 20 counts per bin using grppha. 

In addition to the XMM-Newton spectra, we also obt ained Swift BAT spe ctra. The spectra are 
time-averaged over 22-months and are fully described in iTueller et all (|201ol fl The BAT spectra 



cover the 14-195 keV band and consist of eight energy bins with edges of 14, 20, 24, 35, 50, 75, 100, 
and 150 keV. 



2.2. XMM-Newton Light Curve Analysis 

In order to search for variability in each of the sources, we created background subtracted 
light curves for pn in a soft (0.3-3 keV), medium (3-7 keV), and hard (7-10 keV) band, using a 
bin size of 128 s. The same source and background regions were used as those to extract the 
spectra. For NVSS 193013-1-341047 the average background subtracted count rate in each band 
was 1.88 X 10-2 cts s'^ (soft), 4.79 x lO^^ cts s'^ (medium), and 1.69 x 10"^ cts s'^ (hard). For IRAS 
05218-1212, the average count rate was 5.15 x 10^^ ctss^"*^ (soft), 4.67 x 10^^ ctss^^ (medium), and 
9.98 X 10-3 cts s^i (hard). 

As a simple search for variability, we computed = ~ /"^)/^^ for each light curve, 

where Cj is the background subtracted count rate in each bin, /x is the average count rate, and 
a is the error. For NVSS 193013-^341047, the reduced values are 0.85 (soft), 0.69 (medium), 
and 1.13 (hard), where the light curves had 113 dof. For IRAS 05218-1212, the reduced values 
are 0.77 (soft), 1.09 (medium), and 1.09 (hard), where the light curves had 101 dof. The amount 
of variability is low in all of the light curves. However, the medium band, which encompasses the 
Fe Ka emission region, corresponds to the greatest variability in NVSS 193013-1-341047 and the 
soft band corresponds to the greatest variability in IRAS 05218-1212. The band with the greatest 



^Swift BAT spectra for the 22-month sample are pubHcly available at 
|http : //heasarc . nasa . gov/docs/ swif t/results/bs22mon/ | 
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short-term variability for each source also corresponds to the band with the highest count rate. 
While we do not present the results of the MOS light curve analysis, we find that the MOS results 
are in agreement with the pn. 



2.3. Spectral Fitting 



Spectral fitting was done using XSPEC vl2.4 ([Arnaudl 119961 ) . For each of the sources, we 
simultaneously fit the XMM-Newton EPIC (pn + MOS) spectra with the Swift BAT spectra, 
covering a wavelength range from 0.3-195 keV. As an initial fit to the spectra, we used a power law 
model absorbed by the Galactic column density in the line of sight, in XSPEC tbabs*powerlaw, to 
each source (Nr = 1.73 x 10^^ cm"^ for NVSS 193013+ 341047 and Nr = 9.83 x lO^" cm'^ for IRAS 
05218-1212, as obtained from lDickev &: Lockmanlll99d ). We added a fiux constant to the model to 
account for differences in the flux between the pn spectrum and each of the MOS and BAT spectra. 
This simple model was not a good fit to either spectrum, with /dof = 943/330 (reduced = 
2.9) for NVSS 193013+341047 and 1184/307 (reduced = 3.9) for IRAS 05218-1212. 



In I Winter et al.l (j2008l ) and I Winter et al.l (l2009al ). we describe our choice of a partially covering 
absorption model (pcf abs in XSPEC) to fit more complex spectra in the Swift AGN catalog. This 
model, M{E) = fcover x exp^^HO-(^^) — f^over), assumes that the intrinsic absorbing column 
density only partially covers the AGN emission and that some portion of the direct emission is 
not blocked by the obscuring region. Adding this model to each of our source spectra, we find a 
large statistical improvement in the fit (Ax^ is 127 for NVSS 193013+341047 and 586 for IRAS 
05218-1212). 

In addition to the partia lly covering absorption model, we added a refiection model, pexrav 
(jMagdziarz &: Zdziarskilll995l ). to account for X-ray emission refiected into the line of sight. The 
parameters of this model include the power law index and normalization, which we fixed to the 
values of the direct emission, the refiection factor (R = 0,/2tt, where Q is the solid angle over which 
the material is refiected), the source redshift, cutoff energy, and the inclination angle (which we fix 
at the default of cos 9 = 0.45). Negative values for the refiection factor indicate a pure refiection 

spectrum. Therefore, we allowed this parameter to vary within the range of 5. The addition 

of the refiection model was a significant improvement to the fits for both spectra, with Ax^ of 24 
for NVSS 193013+341047 and 31 for IRAS 05218-1212. However, we could not constrain the cutoff 
energy (instead we fixed this to 100 keV) or the refiection parameters. 

The final component added to our spectral model was a gaussian line (zgauss) to fit the 
Fe Ka emission. A line was statistically significant in both sources, with Ax^ of 17 for NVSS 
193013+341047 and 21 for IRAS 05218-1212. In addition to the Fe Ka emission at 6.41 keV, the 
spectrum of IRAS 05218-1212 is more complex. In Figure [3l we show the ratio of the data/model 
in the 5-8 keV band for this source with the best-fit model, excluding the gaussian added to fit Fe 
Ka. The most significant line is clearly the narrow 6.41 keV emission. Also present is a narrow 
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emission line at 6.0 keV. 



Full details of our best-f it spectral parameters are presented in Table [2j A similar analysis is 
found in iTrippe et alj (|201ll ) for IRAS 05218-1212, which uncovers similar best-fit values as our 
analysis. To summarize our results from the spectral fits, we find that the spectra of both sources 
are heavily absorbed with Nh > 10^'^ cm~^, but Compton thin. The absorption is partially covering, 
covering approximately 93% of the emission in both sources. The measured power law indices of 



both s ources deviate from the average value of 1.75 found for the Swift BAT sample (jWinter et al 



2009al ). with T being fiat for NVSS 193013+341047 (1.42) and steep for IRAS 05218-1212 (2.11). 
Refiection is statistically significant in both sets of spectra, but not well constrained. Additionally, 
an Fe Ka emission line is present in both spectra with an equivalent width of ~ 150 eV. The width 
of this fine is slightly broad, with a = 0.14 keV, for NVSS 193013+341047. 



3. Comparison with Other Bands 



On the basis of the X-ray spectra alone , these two AGN clea rly have complex X-ray spectra. 
With higher signal-to-noise X-ray spectra, IWinter et al.l (|2009al ) found that more than 50% of 
the 153 9-month BAT AGN have complex spectra well-fit by a partial covering model. Among 
these, nearly all were absorbed sources (A^h > 10^^ cm^^). While the initial Swift XRT + BAT fits 
suggested that these sources were Compton-thick candidates, our XMM-Newton observations reveal 
that NVSS 193013+341047 and IRAS 05218-1212 are Compton thin. However, the high column 
densities above 10^^ cm~^ and covering fractions > 90% show that these are likely hidden/buried 
AGN. This is a new class of objects, comprising 24% of the Swift AGN, uncovered in high energy 
surveys, whose AGN nature is hidde n in the optical/soft X-rays due to obscuring gas and dust 
(|Ueda et ahlboOTl : Iwinter et all2009al E]) . The most interesting properties of these sources, however, 
are evident when comparing the X-ray properties with other wavelengths. These include a mis- 
match between the optical and X-ray classifications and unusually red spectral energy distributions 
(SEDs). 



3.1. Optical 



Both the optical spectra of NVSS 193013+341047 and IRAS 05218-1212 show clear broad 
lines, leading to an optical classification of Sy 1-1.8, typically associated with an unabsorbed AGN. 
For NVSS 193013+341047, it is unlikely that the X-ray/optical mis-matc h is due to variab ility 
since the broad line s were evident in two spectra taken over 3 months apart (jLandi et al.l 120071 and 
Winter et al.ll2010al ). Further, Swift XRT spectra, indicating the X-ray spectrum to be absorbed, 
were also taken within a few months of the opt ical. For IRAS 05218-1 212, an EXOSAT spectrum, 
taken within 2 years of the optical spectrum of lMorris WardI (|l989l ). showed the 2-1 keV X-ray 



band to be a factor of 10 times under-luminous compared to other IRAS-selected AGN (jWard et al 



- 6 - 



19881 ). indicating that the source was Ukely heavily absorbed at this time. Additionally, an optical 
spectrum from 22 April 2009, presented in iTrippe et al.l (|201lh . also confirms the optical broad 
emission line spectrum. Thus, the mismatch between X-ray and optical spectra is unlikely due 
to variability in the sou rces. Further, the optical polarization fraction of IRAS 05218-1212 was 
reported as less than 1% (jBrindle et al.lll99d ). showing that the optical light is not mostly scattered 
around the X-ray obscuring region. 

In the XMM large sc ale structure surve y, 7 out of 61 (11%) optical broad line sources exhibited 
X-ray absorbed spectra (jGarcet et al.l 120071 ). It is interesting to note that all of these sources are 
more luminous (in lj2-iokeV by ~ 100) and at m uch higher re d shift (< z >= 1.93) than ou r 
target sources. Similar sources have been found bv iPerola et al.l (j2004l ) and lEckart et al.l (|2006l ). 
comprising 5 - 15% of their X-ray samples. Thus, our targets may be the first low-redshift AGN 
discovered from this class of object, which includes a significant fraction of high redshift AGN. As 
the low-redshift analogs of this class, we will be able to study the sources and hosts, in particular, 
in much greater detail due to their much closer distance. 



3.2. Infrared 



Not only are the X-ray and optical spectra of our targets puzzling, but both sources are very 
red, with J-K = 2.43 (NVSS 193013+341047) and 1.74 (IRAS 05218-1212). These values are high 
compared to the 9-month BAT AGN sample, taken from NED. In Figure HI we plot the 2MASS 
J — K values for the 9-n ionth sample versus both X-ray absorbing column and the 14-195 keV 
luminosity, obtained from lwinter et aP (2009a). These plots show a few interesting results. First, 
Syls (triangles) are redder than Sy2s (circles), with the average J — K of 1.32 for the 31 Sy 1-1. 5s, 
1.13 for the 26 Sy 1.5-1. 9s, and 1.07 for the 30 Sy2s. The high J-K values of NVSS 193013+341047 
and IRAS 05218-1212 align them more closely with the Syls. This also causes our sources, along 
with 3C 105, to inhabit a strange place in the plot oi J — K vs. Nh, being higher column density 
sources with higher J — K values. 3C 105, however, is unlike our sources in that it is a narrow line 
radio galaxy which does not exhibit broad optical lines. 

Another interesting result we find is that the sources with the highest J — K values are also the 
most luminous of the BAT sample. In addition to NVSS 193013+341047 and IRAS 05218-1212, 
3C 105, 3C 111.0, 3C 273, and EXO 055620-3820 are the most luminous sources with the highest 
J — K. Among these, only 3C 105 has a high column density. The host galaxie s of both 3C 105 an d 
3C 111.0 are identified as ellipticals from HST observations with NICMOS (|Madrid et ahlbood ). 
3C 273 is a blazar, while EXO 05 5620-3820 is a Syl with a complex X-ray spectrum best-fit by 
an ionized partial covering model ([Turner et al.l 119961 ). The host types are unknown for both of 
these sources. Both NVSS 193013+341047 and IRAS 05218-1212 are unusual, even among the most 
luminous and highest J — K sources in the sample, in having a mismatch between the X-ray and 
optical types. While the X-ray spectra are similar to 3C 105, this source exhibits no optical broad 
lines. 
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3.3. Broadband SEDs 

In Figure O we show the observed broad-band SEDs for both NVSS 193013+341047 and 
IRAS 05218-1212. The radio/IR data includes publicly available flux densities from NVSS, IRAS, 
and 2MASS (obtained through the NASA/IPAC Extragalactic Database). In addition to this, we 
include our own measurements from Swift BAT and XMM-Newton - including both the observed 
X-ray measurements based on our best-fit model and Optical/UV measurements from the Optical 
Montior (OM) and the Swift UV/optical telescope (UVOT). 

Data from the OM were available in the U and UVWl filters for IRAS 02218-1212 and in the 
UVWl filter for NVSS 193013+341047. We processed the OM data with the SAS script omichain. 
The flux and magnitudes were extracted with the interactive SAS package omsource, using a 
6" circular aperture centered on the source and a circular background region, free of sources. We 
also include Swift UVOT observations for our target sources, which were obtained using the same 
region flies as for the OM data, and extracted with the Swift tool uvot source. Processed Swift 
UVOT observations were available in the V, B, U, UVWl, UVM2, and UVW2 bands for NVSS 
193013+341047 (with Swift observations 00035274002 and 00035274003) and the V, UVM2, and 
UVW2 bands for IRAS 02218-1212 (Swift observations 00037085001, 00037085002, 00037085003, 
and 00037085004). Our broadband SEDs also include the X-ray flux densities measured in the 
0.5-2 keV, 2-10 keV, and 14-195 keV bands, based on our joint XMM + BAT spectral flts. The 
measurements are included in Table [3j 

The spectral energy distributions (SEDs) for both sources, particularly for NVSS 
193013+341047, show the sources to be very bright in the IR with sharp declines in flux into 
the UV and soft X-rays. While bright IR emission could be the result of reprocessing of the AGN 
emission in dust, which obscures the optical/soft X-rays, it is still unclear as to why there are 
optical broad lines. We are unaware of any other low-redshift AGN of type 1 or 2 with similar 
SEDs. While other red AGN {J - K > 2) have been selected through 2MASS, X-ray fol low-ups 



have not shown the same mismatch between X-ray and optical spectra ()Wilkes et al.ll2005l ). 



Discussion 



4.1. Optical and X-ray Reddening 



Multi-wavelength observations of active galaxies previously revealed differences in the redden- 
ing measured between X-ray and optical spectroscopy. Assuming a standard Milky Way extinction 
curve and dust-to-gas ratio, derived X-ray column densit ies can be up to t e n to one hundred times 
that of estimates bas ed on optical emission lines (e.g., iMushotzkvl Il982l : iMaccacaro et al.l Il982l : 
Maiolino et al.l l2001bl ). In the Swift sample, optical estimates of the reddening from the ratio 
of narrow Ha/II/3 emission lines show this same trend with the X-ray derived hydro gen columns 
10-100 times larger than the optically derived column densities (jWinter et al.ll2010bl ). For NVSS 
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193013+341047 in particular, our estim ate of the reddening from the ratio of narrow Ha/H/3 is 
E{B — V) = 0.16 (jWinter et al.ll2010al ). corr esponding to an optic ally derived column density of 
8.5 X lO^'^ cm~^ (based on the conversion from iMaiolino et al.ll2001bl ). Therefore, for this source the 
X-ray derived column is nearly 300 times larger than the optically derived column, on the extreme 
end of what is typically found for differences in the optical/X-ray derived reddening values. 

There are a number of ideas in the literature about how to account for differences in the X- 
ray/optical obscuration. One possible explanation for mismatches in the optical/ X-ray reddening 



i s that the dust-to-gas ratio for AGN is lower than in the ISM in our own Galaxy. IMaiolino et al 



(j2001al ) suggest that large dust grains in the region around the AGN are responsible for the differ- 
ences between t he estimates of the opt ical and X-ray reddening. Among the evidence for unusual 
dust properties, IMaiolino et al.l (j2001bl ) show that active galaxies have faint to no detection of the 
silicate features at 9.7 /xm and the 2175 A carbon feature, both of which indicate the presence of 
small dust grains. Several other possibilities to explain the optical/X-ray mismatch that are dis- 
cussed in these papers include higher metallicities in active galaxies, the hard X-rays being absorbed 
by the optical broad emission line producing clouds (which are expected to have Nh ~ 10^^ cm~^), 
and scattering of the broad line emission t hrough a region of lower column density than what 
obscures the X-ray emission. Alternatively, IShi et al.l (|2006l ) propose an accretion disk geometry 
that can explain the X-ray obscuration and 9.7 //m silicate feature. Yet another explanation for 
differences in reddening is that the X-ray absorption is related to an AGN-driven wind. We describe 
this possibility in the following sub-section. 



4.2. Ionized Emitters 



Mis-matches in classification between the X-ray absorbed and broad-line optical spectra are 
sometimes explained by the presence of a high column of ionized gas. Such a model has been adopted 
to account for active galaxies at both low and high redshift. For instance, the optical and UV spectra 
of the low-z AGN NGC 4151 reveal broad emission lines while the X-ray grating observations 
show that the soft X-rays are heavily obscured and dominated by strong ionized emission (e.g., 
Kraemer et al.ll2005l ). NGC 4151 is classified similarly as a Sy 1.5, with a strongly de tected multiple- 
component ionized outflowing wind in the UV/X-rays (jKraemer et al.ll2005l . l2006l ). NGC 4151 is 
also detected in the Swift catalog, however, both the amount of reddening (from J-K) and the 
column density of gas are much lower than the two new sources we present in this current paper, 
showing that it is not a useful analog for heavily reddened NVSS 193013+341047 and IRAS 05218- 
1212. 



Perhaps a more relevant example of a source with a mis-match between optical and X-ray 
spectra is the broad-line radio galaxy 3C 445, which is at a similar redshift to our sources with z = 
0.057. Broad emission lines are observe d in this source (FWH M ~ 6400 km s~^) along with a highly 
reddened continuum (Eb-v = 1 mag; ICrenshaw et al.lll988l ). Suzaku and Chandra observations 



show that 3C 445 is heavily absorbed with a column density of ~ 10^^ cm and that the soft 
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spectr um, like that of NGC 4151, is dominated by emission lines (jBraito et al.ll201ll ). iBraito et al 
(|201ll ) conclude that the photo- ionized emitting gas is likely associated with the BLR gas, which 
we are able to view because the high column density absorber is not completely blocking the BLR. 
This is possibl e if the obscuring r egioi i is not a uni form torus but instead clumpy or wind- like in 
structure fe.g.. iRisaliti et alJbooa and lElitzuilbood ). 



In both NVSS 193013+341047 and IRAS 05218-1212 the residuals in the soft X-rays of our 
spectra suggest that emission lines may be present, but our observations do not have the resolution 
to confirm them. If higher resolution X-ray observations indicate both the high column density 
neutral obscuration and a strong emission component, then the same clumpy torus/wind models 
used to explain NGC 4151 and 3C 445 also explain the mis-match in classifications for our unusual 
sources. To confirm strong emission, potentially associated with an outfiow, we will propose for 
grating observations of these sources with XMM-Newton/ Chandra. 

Recently, active galaxy/quasar driven ou tflows were also invoked to explain highly absorbed 
QSOs detected with ROSAT/sub-mm surveys (jPage et allboilh . In this case, the z ~ 2 QSOs with 
rest-frame broad UV emission lin es and heavily absorbed X-rays have X-ray spectra that can be 
modeled with ionized absorption. [Alexander et al.l (|2005l ) suggest that these sources may be in a 
transition between the absorbed and unabsorbed phases. As lPage et al.l (|201ll ) show, this stage may 
be marked by strong AGN-driven winds. If this is the case, then the intermediate Seyfert 1 types, 
like the newly detected Swift sources we discuss in this paper, are likely in a similar transitional 
stage. This is supported by the fact that many of the well-known Seyferts with strong outfiow 
detections are Sy 1.5s (e.g., NGC 3516, NGC 3783, and NGC 4151). 



4.3. Low-redshift EROs? 



While there are no known low-redshift AGN with similar SEDs to NVSS 193013+341047 
and IRAS 05218-1212, the high redshift universe does provide a useful analog. Namely, we find 
that the SEDs of our sources, particularly NVSS 193013+341047, match that of XBS J0216-0435 
^ 2). This source is identified as a type 2 quasar hosted in an ERG, whose SED is presented in 
Severgnini et al.l (j2006l ). In Figure El we include the SEP of thi s source with a cross symbol. As 
pointed out and illustrated in Figure 6 of Severgnini et al. (j2006l ). the SED is unlike those of radio 
quiet or radio loud AGN. 



EROs are identified as red sources with R — K > 5 (jElston et al.lll988l ) oi I — K > 4. Joint 
IR/X-ray st udies of EROs have fou nd that most AGN in this category are heavily obscured (A^h > 



lO^^cm"^) (lAlexander et al.ll2002l ). Based on the 2MASS (J, H, K) colors, we estimate R - K 



of 4.8 (NVSS 193013+341047) and 3.5 (IRAS 05218-1212). This places them below the formal 
definition of an ERO. However, with a similar SED, very red colors, and X-ray absorbed spectra, 
our sources may be the low redshift equivalent. Particularly, the result showing J — K values to 
increase with increasing luminosity (Figured]), supports this claim, as our sources are a factor of 
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10-100 less luminous in the 2-10 keV X-ray band. If these sources are low redshift EROs, then 
multi-wavelength studies of their properties will help to answer many of the questions surrounding 
these high redshift sources. 

Previous X-ray observations of EROs sugges t that they are absorbed sources with Nh > 



10 cm , representing a class of type II quasars (jMignoli et al.ll2004l : iBrusa et al.ll2005al ). These 



type II quasars all tend to be hosted within elliptical galaxies. While we can not be sure of the 
morphology of NVSS 193013+341047 and IRAS 05218-1212, no spiral arms are distinguished in 
the available Digital Sky Survey images. The morphologies of Swift-detected AGN with similar 
J — K colors all tend to be in ellipticals, so it is possible that our unusual sources are also hosted 
in ellipticals. If so, our highly reddened active galaxies are two of only a handful of Swift-detected 
AGN with el liptical hosts, as the majority of the very hard X-ray detected galaxies are associated 



with spirals (IKoss et al. 



2011 



5. Summary 

In this paper, we present two very unusual AGN discovered by the Swift BAT all-sky survey. In 
the X-ray band, these sources are clearly heavily absorbed (Nh =1-2 xlO^^ cm~^), with little X-ray 
flux below 2 keV. Both of these sources exhibit broad lines in their optical spectra, seemingly at odds 
with the heavily absorbed X-ray spectra. Further, the SEDs, revealed through publicly available 
NVSS, IRAS, 2MASS, and Swift data, are very red. As such, these sources appear to be the low- 
redshift analogs of AGN hosted in EROs. It is still unclear, however, how the multi-wavelength 
properties fit together, particularly the mismatch between X-ray and optical classification, but one 
possible explanation is that the heavy obscuring column density is only partially obscuring the 
BLR, as well as an associated X-ray emitting region associated with an outflow. High signal-to- 
noise grating observations are required to conflrm outflowing ionized gas in the spectra of these 
sources. 

Another possible explanation for the mis-match in classiflcations is that the dust properties in 
the host galaxies of NVSS 193013+341047 and IRAS 05218-1212 are different. Infrared spectroscopy 
could help to uncover the dust properties in these sources. What we can discern, on a larger scale, 
is that the host galaxy types of these unusual sources may be different from the majority of the 
BAT-selected sample. From the 2MASS derived absolute magnitudes of J = —22.8 and —23.78, 
we flnd that the host galaxies are of typical luminosity (L*). Unfortunately, the available images 
from UVOT and the Digital Sky Survey, are not of sufficient quality to resolve structure, determine 
whether the sources are actively undergoing star formation, or discern dust lanes. Since spiral arms 
are not indicated and two of the BAT AGN with high 14-195 keV luminosities and similar J — K 
values are hosted in ellipticals, it is possible that the hosts of these sources are also ellipticals. If 
this is the case, the sources are particularly interesting, being among only ~ 5/100 BAT AGNs 
known to be hosted in an elliptical. 
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Finally, and possibly most importantly, we showed that the 2MASS-derived J — K values for 
the BAT AGN appear to be connected to both the amount of X-ray obscuration and the 14-195 keV 
luminosity. We showed that J —K \s lower for Sy2s and higher for Syls as well as for high luminosity 
sources. This is interesting considering that at high redshift { z ga 1-2) the redde st EROs are now 



being used as a means to select obscured Type II QSOs (e.g., iBrusa et al.ll2005bl ) 
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Table 1. Details of the XMM-Newton Observations 



Soiirce 


RA 


Dec 


z 


Obs. ID 


Date 


Exp. Timet 


<Ct Rate>t 


NVSS 193013+341047 


19:30:13.86 


+34:10:50.04 


0.0629 


0602840101 


16-05-2009 


12.91 


0.29 


IRAS 05218-1212 


05:24:06.80 


-12:10:10.87 


0.0490 


0551950401 


24-08-2008 


11.27 


0.17 



tThe exposure time (in units of ks) and count rate (in units of ctss ^) are recorded for the pn spectra in the 0.1-10 keV 
band. 
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Table 2. Joint XMM-Newton and BAT spectral fits 



Parameter NVSS 193013+341047 IRAS 05218-1212 





24.5tli 




J cover 


0.9310:°^ 


0.93tni 


rb 


1.42l°:J° 


9 1 1 +0.09 
^•■'■■'■-0.10 


At" 


9.85 xlO-^ 


1.74 xlO-3 




(-5) 


(-5) 




6.381°;°^ 


6.4112;^^ 


a'' 




0.07t«;°? 




149t^^ 


1 79 + 137 
'^'■^-63 


xVdof 


292/323 


302/300 



'^Parameters of the pcf abs model, including hiydrogen col- 
umn density measured in units of 10^^ atoms cm~^ and the 
covering fraction of the absorber. 

^'Spectral index computed with the reflection model pexrav 
assuming total reflection and a fixed energy cutoff at 150 keV. 
The reflection parameter is not well-constrained. The quoted 
normalization on the power-law component is in units of pho- 
tons keV~^ cm~^ s~^ at IkeV. 

'^Equivalent width and normalization (with errors or as an 
upper limit) on a fixed Gaussian Fe K line at 6.41 keV in units 
of eV and 10~^x total photons cm~^ s~^, respectively. 



-15- 



Table 3. Multi- wavelength Flux Measurements 







iNVoo iyoUio+o4iU4 / 


Tt? A Q nP;01 Q 1010 


Measurement 


u [tizj 


(.mJyJ 


(Jyj 


NVSS 


1.4x10^ 


4.4 ±0.6 


9.2 ± 1.1 


IRAS 100/x 


3x10^^ 




1020 


IRAS 60 fj. 


5x10^^ 




456.3 ± 50 


IRAS 25/Li 


1.2x10^^ 




280 


IRAS 12 ij, 


2.5x10^^ 




100 


2MASS K, 


1.388571x10^'^ 


12.6 


16.5 


2MASS H 


1.803805x10" 


6.05 


11.3 


2MASS J 


2.427469x10" 


3.2 


7.95 


UVOT B 


6.662055x10" 


0.277 ±0.049 




UVOT B 


6.662055x10" 


0.26 ±0.01 




OM U 


8.30781414x10" 




2.41 ±0.10 


UVOT V 


8.446x10" 


0.78 ±0.09 


6.15 ±0.10 


UVOT U 


8.81742524x10" 


0.13 ±0.02 




UVOT U 


8.81742524x10" 


0.196 ±0.019 




OM UVWl 


9.82092118x10" 


0.05 ±0.01 


2.10±0.01 


UVOT UVM2 


1.297803x10^5 


0.027 ±0.006 


1.49 ±0.01 


UVOT UVM2 


1.297803x10^5 


0.024 ±0.007 


1.64 ±0.02 


UVOT UVW2 


1.414115x10^5 


0.030 ±0.005 


1.41 ±0.02 


UVOT UVW2 


1.414115x10^5 




1.54 ±0.02 


PN soft 


1.25643545x10^^ 


0.00010 


0.00036 


PN hard 


1.89640880x10^8 


0.00015 


0.00027 


BAT 


6.01787337x10^8 


0.00013 


0.00029 



Notes: The NVSS, IRAS, and 2MASS flux densities were obtained using the NASA/IPAC Extragalactic 
Database. Swift UVOT and XMM-Newton OM observations were analyzed as explained in the text. The 
X-ray measurements, obtained from the best-fit X-ray models, correspond to the 0.5-2 keV soft band, 2-10 

keV hard band, and the 14-195 keV BAT band fluxes. 
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Fig. 1.— The XMM-Newton pn and Swift BAT spectra are plotted for NVSS 193013+341047, 
along with the best fit model, as well as the ratio of the data/model in the bottom plot. The 
spectrum is complex and required a partial covering absorber model. 
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Fig. 2.— The XMM-Newton pn and Swift BAT spectra are plotted for IRAS 05218-1212, afong 
with the best fit model, as well as the ratio of the data/model in the bottom plot. The spectrum 
is complex and required a partial covering absorber model. 
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o ' ' ' ' 

5 

Energy (keV) 

Fig. 3. — The ratio of the data/model for a partially covering power law model without accounting 
for the Fe Ka emission is shown in the 5-8 keV band. The pn (black), MOSl (red), and M0S2 
(green) data points are shown. In addition to the narrow Fe Ka line at 6.41 keV, a second line is 
significant in each of the detectors (Ax^ = 9) at 6.0 keV. 
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Fig. 4. — Plot of reddenning versus the X-ray measured column density (Nh) and the 14-195 keV lu- 
minosity from BAT. The reddenning (J-K) is constructed from the 2MASS values, publicly available 
in NED. NVSS 193013-^341047 and IRAS 05218-1212 (pentagons) are among the reddest sources, 
along with 3C 111.0, EXO 055620-3820, 3C 105, and 3C 273. Seyfert 1-1.5 sources (triangles), 
Seyfert 1.5-2 (squares), and Seyfert 2s (circles) are all indicated. 
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Fig. 5. — We constructed the observed spectral energy distributions (SEDs) for our targets using 
available radio/IRAS/2MASS measurements as well as our own measurements from the SWIFT 
and XMM-Newton observatories (including optical/UV from UVOT and XMM OM, 0.5 - lOkeV 
X-rays with XRT/XMM pn, and 14 - 195 keV very hard X-rays with BAT). Both sources show a 
steep decline in flux from infrared into ultraviolet and soft X-rays, not typical of the other BAT 
AGNs. However, the sharp d rop in flux from the IR to UV/soft X-rays is similar to the SED of 
ERO source XBS J0216-0435 (|Severgnini et al.l 120061 ) . plotted with the cross symbols. 



